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ABSTRACT 



Fast coronal mass ejections (CMEs) generate standing or bow shocks as 
they propagate through the corona and solar wind. Although CME shocks have 
previously been detected indirectly via their emission at radio frequencies, direct 
imaging has remained elusive due to their low contrast at optical wavelengths. 
Here we report the first images of a CME-driven shock as it propagates through 
interplanetary space from 8R Q to 120 R (0.5 AU), using observations from 
the STEREO Heliospheric Imager (HI). The CME was measured to have a 
velocity of ~1000kms _1 and a Mach number of 4.1±1.2, while the shock 
front stand-off distance (A) was found to increase linearly to ~20 R Q at 
0.5 AU. The normalised standoff distance (A/ Do) showed reasonable agreement 
with semi-empirical relations, where Dq is the CME radius. However, when 
normalised using the radius of curvature, A/Rq did not agree well with theory, 
implying that Ro was under-estimated by a factor of i=s3-8. This is most likely 
due to the difficulty in estimating the larger radius of curvature along the CME 
axis from the observations, which provide only a cross-sectional view of the CME. 



Subject headings: Sun: coronal mass ejections (CMEs) — Shock waves 



1. Introduction 

Bow shocks occu r when a blunt object moves relative to a medium at supersonic speeds 
( iRathakrishnanl 120101 ) . These shocks are formed across many sc ales and in different co ndi- 
tions; from astrophysical shocks such as p lanetary bow shocks (Slavin fc Holzerl 1 198 if ), or 
the shock at the edge of the Heliosphere (Ivan Buren et al.lll995l ). to shocks generated by 
the reentry of the Apollo mission capsules ( Glass! 1977 ). Coronal mass ejections (CMEs) 
which travel faster than the local fast magnetosonic velocity (with r espect to the solar w ind 
velocity) produce such standing shocks in the frame of the CME ( IStewart et al.lll974al Jbl). 
Interplanetary (IP) CME-driven shocks have previously been detected in radio observations 
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as Type II bursts and using in- situ measurements. Direct im aging of shocks, on the other 



hand, has remained elus ive, primarily due their low contrast ( IVourlidas Ontiverosl 12009 
Gopalswamy et al.ll2008l ). 




Fig. 1. — Diagram of the various quantities used to describe the shock and CME. 

The shape, size, and standoff distance of a shock are controlled by several factors: the 
shape and size of the obstacle; the velocity difference between the obstacle and the medium 
with respect to the sonic speed (i.e., the Mach number); and the properties of the medium, 
such as the ratio of specific heats (7) and the magnetic field. Relationships between the 
shock standoff distance and the Mach number have be en derived by a number of different 
authors. The well known semi-empirical relationship of ISeiffl ( 119621 ) has the form 

A = 0.78^, (1) 
Do Pd 

which was derived for a spherical object, where A is the shock standoff distance, Dq is 
distance from the centre to the nose of the obstacle, in this case the radius, and p u , pd are 
the densities upstrea m and downstream of the shock respectively. Using gas-dynamic theory, 
Spreiter et al.l (119661 ) demonstrated that p u /pd could be written in terms of the the upstream 
sonic Mach number, M s , and the ratio of specific heats 7: 



A 



1.1 



( 7 -l)M s 2 + 2 



(2) 



( 7 + l)M2 • 

The increase in the coefficient in the standoff relations from 0.78 to 1.1 is due to the fact 
the object under consideration (Earth's magnetosphere) in Equation (J2J) is more blunt than 
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a sphere; specifically, it is an elongated ellipse. Neither Equations (00) or fl2]) behave as 
expected at low Mach numbers, where the shock should move to a large standoff distance. 
A modification which corrects for this enables Equation ([2]) to be written in the form 



A 



(7 - 1)M 2 + 2 



• (7+l)(M 2 -l)' w 

where the additional term in the den ominator ensures the s hock moves to a large distance as 
the Mach number approaches unity ( iFarris fc Russelllll994j ). They also suggested that using 
the obstacle radius of curvature rather than radius would be more suitable as it accounts for 
the shape of the obstacle, resulting in 



A ( 7 -l)M s 2 + 2 

Ro ' (7+l)M 



1) 



(4) 



where i?o is the obstacle radius of curvature. 



In general, a conic section can be represented by y(x) 2 = 2R(D — x) + b(D — x) 2 where 
b is the bluntness (6 < —1: blunt elliptic; b = —1: spherical; — 1 < b < 0: elongated 
ecliptic; b = 0: parabolic; 6 > 0: hyperbolic). The shape of the shock fronts are known to 
be re presented by a mod ified conic section. One such parameterisation of the shock front, 
from IVerigin et al.l ( 120031 ). is: 



y 2 (x) 



2R S (D 



x) + 



x) 



Ml - 1 
6 5 (M 2 -1)-1 , 



v ^ 1 l + d s {D s -x)/R s " (5) 

where 65 is the bluntness of the shock, and d s is related to the asymptotic downstream slope 
or Mach cone (see Figure [T|) . 

The relationships between the standoff distance and the Mach number have been investi- 
gated from a number of perspectives, including numerical modeling, anal ytical relations, lab- 
oratory exp eriments and in- situ measurements of planetary bow shocks (jSpreiter &: Stahara 
19951 . Il980l ). These have shown that in general the semi-empirical relations p rovide an ad- 



equat e description of shocks, with the low Mach regime being an exception (jVerigin et al. 



20031 ). Depending on the physical, context the sonic Mach number (Ms) can be replaced with 
the magnetosonic Mach number (Mms), when dealing with plasmas such as the solar wind 
and CMEs. It has been shown that using gasdynamic relations when dealing magnetised 
plasmas works well when the MHD mach numbers are high. It also provides a good ap- 
proximation when the Alfven (Ma) or fast magnetosonic (Mms) Mach numbers are low and 



these Mach numbers are substituted for the gasdynamic (Ms) Mach numbers (IFairfield et al. 
200lh . 
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Fig. 2. — Simultaneous STEREO observations of the CME and shock front from Ahead (top 
row) and Behind (bottom row). The CME and shock are indicated on the individual images 
where applicable. No CME or shock is visible in the HI 1 A observation (the Sun is off the 
right hand edge of the HI 1 A image). These images have been severely clipped and smoothed 
to make the shock more discernible. The shock is clearly visible in the accompanying movies 
(Ml, M2). 

Standoff distances of CM E-driven shocks haye been investigated from an in- situ per- 



spect iv e by many authors ( e .g.. IRussell &: Mulliganll2002l ; iLepping et al.ll2008l ; lOdstrcil et al. 



2005 ). Russell &: Mulligan ( 20021 ) found the sh ock standoff di s tance (A; thickness of mag- 



netosheath) was of the order of 21 R Q at 1 AU. ILepping et all (120081 ) derived an average A 
of about 8 R at 1 AU. However, when considering the CME radi us (flux rope radius) as Dp 
the typical A expected from Equation (3) is about 5 R at 1 AU. IRussell fc Mulligan! (120021 ) 
proposed that Equation (T4J) may be more suited as it accounts for the fact the CME front 
may not be circular and that the radius of curvature at the nose is a dominant factor in 
determining the standoff distance. However, they found that Equation (p|J did not fit the 
observations either and speculated this may be due to observational effect of only measuring 
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one of the radii of curvature of the CME. The underlying structure of a CME is believed 
to be a flux rope, which has two characteristic curvatures; a smaller one due the curvature 
perpendicular to its axis (the radius when viewed as a cross-section) and the larger curvature 
along the axis. 

In this Letter, we investigate if the shock relations above hold for a CME-driven IP 
shock. Specifically, we use direct observations of a CME-driven shock observed in C0R2 
and HI1 instru ments of the Sun Ear th Connection Coronal and Helios pheric Investigation 
suite (SECCHI, iHoward et all l2008h on STEREO kaiser et al. l l20Q8h . In Section EJ we 
present SECCHI observations of the CME and resulting shock and describe the analysis 
technique. The results of our analysis are presented in Section [3) We discuss our results and 
state our conclusions in Section @] 



2. Observations and Data Analysis 



The CME analysed first appeared in the COR1 ([Thompson et al.l 120031 ) coronagraph 
images from STEREO B (Behind) at 15:55 UT on 2008 April 5. It was most likely associated 
with a B-class flare from NOAA active region 10987, which was just behind the west limb 
as viewed from Earth. Figure [2] shows the CME as it propagates out from the Sun into the 
different instruments' field-of-views from 8R to 120 R Q . The CME was visible in both A 
(Ahead) and B spacecraft in the i nner and outer cor onagraphs (COR1 and COR2; see movie 
Ml), but was only visible in HI1 fevles et al-lbopgh from STEREO B (see movie M2). The 
CME propagation direction was found to be ~106° west of the Sun-Earth line, the spacecraft 
were at a separation angle of 48° degrees (from each other). The shock is visible as a curved 
brightness enhancement in both the COR2 images in Figure [2] and also the HI1 B image. 
The accompanying movies Ml and M2 show the shock more clearly. We have made the 
assumption that the curved front is a shock, there are no radio or in- situ data available to 
corroborate this. However, due to the CME's velocity (~1000kms -1 ) and the smoothness 
and position of the feature ahead of the CME, it can be argued that this it is a legitimate 
assumption (IBemporad fc Mancusoll2010l ; lOntiveros fc Vourlidasl 120091 ). 



The observations were reduced using secchLprep from the SolarSoft library (IFreeland fc Handy 



19981 ). This corrects for a number of effects such as bias, flat-field and distortions. The coro- 
nagraphs take sequences of three polarised observations, which are combined to produce 
total brightness images. The COR1 and COR2 observations were used to produce standard 
running difference images. The HI observations were background subtracted. Modified run- 
ning diffe rence images were th en created, which account for the motion of the back-ground 
star field jMalonev et alJboOflh . 
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Fig. 3. — (a) 3D reconstruction of the CME and shock front viewed perpendicular to the 
propagation direction with the initial ellipse fits over plotted, (b) Data transformed into a 
coordinate system centred on the initial ellipse fit. Over plotted is the subsequent fit to the 
shock front using Equation (jSJ). The data have been normalised with respect to Do- 



For each observation in which the CME or shock was visible the front was identified. A 
number of points along this front were then manually chosen. For the observations where 
the CME or shock was observed from both spacecraft, the front was localised in three 



dimensions using th e tie-point method (jlnhesterll2006l ; iByrne et al.ll2010l ; iMierla et al.ll2010 



Temmer et al.l 120091 ) . As the CME or shock was only observed by one spacecraft in the 
HI field-of-view, we used the additional assumption of pseudo-radial propagation, based o n 
the direction derived from COR1 and COR2 to localise the front (IMaloney et al.l 120091 ). 
The resulting data consisted of a series of points in 3D for the CME and shock for each 
observation time. Figure |3]Ja) shows the 3D reconstruction of both the shock front and CME 
front, viewed perpendicular to the direction of propagation (assumed to be a cross-section). 
The techniques for deriving the 3D coordinates of features in the COR1/2 and especially 
the HI field-of-view are not without error. In the case of the event studied here, the CME 
was close (<10°) to the plane-of-sky of STEREO A. As a result, errors in position should be 
small. 

In order to compare with relations in Section [TJ the data was transformed into a co- 
ordinate system centred on the CME. To accomplish this, each CME front was fit with an 
ellipse. The centre coordinates of these fits were then used to collapse all the data on to a 
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common coordinate system centred on the CME. The shock front was fit with Equation (JSJ), 
which gave the shock properties such as the shock standoff distance A, the Mach number M, 
and the radius of curvature at the nose of the obstacle Ro- Figure E^a) shows data and the 
initial fit, Figure |3^b) shows the shifted data and the shock fit using Equation fl5J). The fast 
magnetosonic Mach number was calculated using M ms = (v cme — v sw )/v ms , where v cme is the 
CME velocity, v sw is the solar wind velocity and v ms is the fast magnetosonic speed. Since 
v sw and v ms were not known at the position of the CME, a model corona was used to evaluate 
them. This was based on the Parker solar wind solution with a simple di polar magnetic fiel d 
of the form B(r) = Bo(R & /r) 3 , where Bq was 2.2 G at the solar surface ( IMann et al.ll2003l ). 
For each of the paired CME and shock observations the standoff distances A (=D$ — Do) 
were obtained by three different means: (i) using the 3D coordinates of the furthest point 
(max(h), where h = \J x 2 + y 2 + z 2 ) on the shock and the CME as h s hk and h cme respec- 
tively, (ii) the previous method can be applied but to the data in the common coordinate 
system which gave D Q and D s , and (iii) the front fitting procedure also produced standoff 
distances. However the results of method (i) cannot be used with the relations from Section 
[I] as they are not in a CME/obstacle centred coordinate system, but the results from method 
(ii) and (iii) can be compared to Equations (J2J), (G2) and (BJ. 



3. Results 



A summary of the shock properties derived from the observations is shown in FigureH^a)- 
(f) as a function of time. With the exception of the CME (h cme ) and shock heights (h s hk), 
all the properties have been derived from the data collapsed on to a common coordinate 
system with respect to the CME. The gap between the first three data points and others is 
a result of both the CME and shock leaving the COR2 field-of-view and entering the HI1 
field-of-view. The contrast between shock and background in the first three and last three 
observation is extremely low, making identification of the shock difficult. As a result, these 
points are not reliable, and should be neglected. Figure H(a) shows the derived heights of 
the CME and shock as they were tracked from 8 R Q to 120 R (0.5 AU). Using a linear fit to 
hshk — h cme (=A) versus h cme (not shown), the extrapolated standoff distance at Earth was 
found to be ~4OR . Figure IU(b) shows the distance to nose of the CME (Do) and shock 
(Ds) front by (filled symbols), also shown are the values derived from fits to the shock and 
CME front (hollow symbols). The increasing offset between the two is due to their differing 
centres of the coordinate systems, as one is elliptic and the other is parabolic. Figure 0Jc) 
shows the standoff distance A derived using Do and Ds (filled symbols) and from the fits 
to the fronts (hollow symbols). Both are in general agreement and show an increase with 
time. The standoff distance normalised using Do is shown in Figure Hfd). The normalised 
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standoff distance is roughly constant with a mean value of 0.37 ±0.09. The standoff distance 
normalised to the radius of curvature at the nose of the CME (Ro) is shown in Figure IHe). 
The curvature could only be derived from the front fitting, as such only hollow data points 
are shown. Figure H]T ) then shows the magnetosonic Mach number (Mms) derived using: 
(i) the CME speed in conjunction with the coronal model (filled symbols) and (ii) the shock 
front fitted using Equation 5 (hollow symbols). The mean Mach number from the coronal 
model was 3.8±0.6, while a value of 4.4±1.6 was found using the front fitting method. The 
mean Mach number from both methods was 4.1±1.2. 

Figure 5(a) shows the relationship between the normalised standoff distance (A/ Do) 
and Mach number (Mms) f° r a number of models. The Mach numbers were calculated using 
the coronal model (filled symbols) and front fitting (open symbols). The normalised standoff 
distances were calculated using measured values of Dq and D s (filled symbols) and fits to the 
CME and shock fronts (hollow symbols). Both s how good ge neral agreement between our 



observations and the models (<20%). The model of ISeifq (119621 ) shows the poorest agreement, 
although this is not unexpected as it was derived for a circular obstacle and the CME is quite 
blunt compared to a circle. Figure 5(b) shows the relationship between the standoff distance 
normalised by the radius of curvature of the CME (A/Rq) and Mach number for a number 
of models. In this case, Ro can only be derived from the front fitting. These values are 
then plotted as a function of the Mach numbers derived using both methods described above 
(hence, each value of A/D Q appears twice). Our results do not agree with the expected 
relation (Equation (|3J)) and indicate the the radius of curvature Ro is underestimate by a 
factor of ~3-8. One possible reason for this is that we have not considered the effect of the 
magnetic field of the CME and solar wind effects on the shock. However, one would expect 
if this had a significant effect it would also affect the other relation. It should be noted that 
the fast magnetosonic velocity and sonic velocity calculated from our model differ by less 
than 7% after excluding the first three data points as mentioned earlier. This also suggests 
that the magnetic field should not play a maj or role. A more like l y rea son is due to an 



observational affect similar to that suggested by lRussell &: Mulligan! (120021 ). where only one 
radius of curvature of the CME is observed. The observations provide a cross-sectional view 
of the CME along one of its axes. As a result, we have no information on the curvature 
along other CME axes. 



4. Discussion and Conclusions 



For the first time, we have imaged a CME-driven shock in white light at large distances 
from the Sun. The shock was tracked from 8 R Q to 120 R Q (0.5 AU) before it became too faint 
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to be identify unambiguously. The CME was measured to have a velocity of ~1000kms~ 1 
and a Mach number of 4.1±1.2, while the shock front stand-off distance (A) was found to 
increase linearly to ~20 R Q at 0.5 AU. The normalised standoff distance (A/ Do) was found 
to be roughly constant with a mean of 0.37±0.09 . The normalised standoff distance derived 
using Do and D$ and its relation to the Mach number (Mms) were compared to previous 
relations and showed reasonable agreement. The normalised standoff distance (A/ Do) and 
Mach number were also derived by fitting the CME and shock front, which agreed well with 
theory and our other method of estimation. The fitting also allowed us to find the CME 
radius of curvature (Ro) enabling us to test the relationship between A/Rc and the Mach 
number. In this case, the derived ratios did not agree with the theoretical predictions and 
showed a significant deviation. 

The faint nature of the shock front made its identification challenging, and thus, the 
front location and characterisation showed some scatter (Figure H]). For example, the Mach 
numbers in Figure H^f ) show a large amount of variability especially from the front fitting. 
The standoff distances in Figure H^c) show the same trend and the two different methods give 
similar results. In should be noted that the first three and last three data points show a large 
deviations from the rest of the data for a number of derived properties. These correspond 
to very low contrast observations, and hence should be ignored. The Mach number derived 
from our coronal model and CME position and speed, and from the shock front fitting both 
agree. This is a good indication that our methods accurately describe the shock even in the 
presence of large uncertainties. 

Both sets of data for the normalised shock standoff distance A/ Do versus Mach number 
(Mms) derived directly and from front fitting show good general agreement (Figure |5]^a)). 
The standoff distance normalised by the CME radius of curvature (A/Rq) verses Mach 
number (Mms) from either the fits or derived directly do not agree with any of the relations 
(Figure E3^b)). Assuming that a CME can be modeled as a flux rope, it should have two 
radius of curvatures. Our observations are a measure of a combination of these, which 
depends on the orientation of the flux rope. This observational affect implies we may only 
be measuring the smaller of the two and this leads to the underestimation of Ro- Finally, 
the general agreement between the Mach number derived from our model and the Mach 
number derived from the fits, suggest the fitting is not the source of the problem. Using 
a mean Mach number of 4 to get a value of 0.26 for A/Ro ratio and the standoff distance 
calculated at Earth (40 R©), we can estimate the radius of curvature of this CME at Earth 
to be 150 R Q (0.7 AU). 

Imaging observations of CME-driven shocks opens up a new avenue for studying their 
fundamental properties. This type of observation will be highly complementary to radio 
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and in- situ measurements. A complete picture of the shock could then be constructed and 
the derived properties from the different observations could be compared and contrasted. 
Furthermore, the analysis presented here will be applicable to future observations of shocks. 
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Fig. 4. — Shock properties derived directly from the observations (filled symbols) and from 
fits to the shock and CME (hollow symbols) as a function of time, (a) The maximum height 
of the CME front (triangles) and shock front (circles), (b) The distance to front of CME 
(Do) and shock (Ds) in CME centred coordinate system, (c) The shock standoff distance 
A. (d) The normalised standoff distance (A/ Do)- (e) Standoff distance (A) normalised by 
the radius of curvature of the CME (Ro)- (f) The Mach number (M) derived the CME 
velocity and model for the corona (filled circles) and from the fits to the shock front (hollow 
circles). 
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Fig. 5. — (a) Shock standoff distance normalised to Do as a function of Mach number, (b) 
Shock standoff normalised by Ro as a function of Mach number. Also shown are the results 
of a number of semi-empirical models. Filled symbols indicate values derived using a coronal 
model, while hollow symbols indicate values derived using front fitting. 



